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We investigate the Gross-Neveu model on the lattice at finite temperature and chemical potential 
in the limit of an infinite number of fermion flavours. We check the universality of the continuum 
limit of staggered and overlap fermions at finite temperature and chemical potential. We show 
that at finite density a recently discovered phase of cold baryonic matter emerges as a baryon 
crystal from a spatially inhomogeneous fermion condensate. However, we also demonstrate that 
on the lattice, this new phase disappears at large coupling or in small volumes. Furthermore, we 
investigate unusual finite size effects that appear at finite chemical potential. Finally, we speculate 
on the implications of our findings for QCD. 
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1. Introduction 

In order to understand the phases of matter at finite temperature and density it is necessary to 
understand the properties of the non-perturbative vacuum of QCD or related models. The Gross- 
Neveu model ||l|] resembles QCD in many respects and can be solved analytically in the limit of 
an infinite number of fermion flavours N. Here we report on our investigation of this model on a 
discrete space-time lattice at finite temperature and density in the large-Af limit. 

2. The Gross-Neveu model 

Let us start with writing the Euclidean Lagrangian density of the Gross-Neveu (GN) model [[[]] 
in 1 + 1 dimensions, 

JSf = f y>«(x)^«(x)-^ ( £ fWfw) , (2-1) 

a=l L \a=l / 

where \\f a (x) are 2-component Dirac spinors and a is a flavour index. Usually, one introduces a 
scalar field o(x) conjugate to Yla=i Xff a (x)ljf a (x) in order to transform away the 4-fermion coupling, 

Jgf = £ f(i)^ a (i) + i(j(i) 2 + (7(i) f f a (x) ¥ a (x). (2.2) 
a=l # a=l 

The GN model shares many interesting properties with QCD. In particular, it is renormalisable 
and asymptotically free, with the j8 -function given to lowest order by j3(g) = — ^2x~g 3 + 0(g 5 ). 
Moreover it enjoys a 0{2N) x Y global symmetry, where Y is the discrete chiral symmetry 

r : Y^YsY, Y^-YYs, o->-a, (2.3) 

and it exhibits spontaneous breaking of this discrete chiral symmetry which in turn leads to the 
fermions acquiring a non-vanishing mass Oq = (o) (dimensional transmutation) 1 . 

In the large-./V limit where the number of fermion flavours N is taken to infinity while keeping 
A = g 2 N fixed, the model can be solved analytically. One can integrate out the fermions to obtain 
Z = J [da] exp {— 5 cff } with 

S ai = Nljdx^--Trlog[^ + a]\. (2.4) 
The minimum of the effective potential is given by a set of equations, 

do(x)S«/N = ^ " <Vv)Trlog [# + a] = 0, Vx, (2.5) 
and for a homogeneous condensate a(x) = a this set reduces to a single equation 

£ = a ff Trlog0 + (7], (2.6) 



Note that there is no Goldstone boson since T is a discrete symmetry. 



2 



New baryon matter in the GN model 



Urs Wenger 







£7 = 


0.4 










ft 


0.2 




; crystal phase 


n 







0.4 



0.8 

77 M/Ob 



1.2 




Figure 1: Revised phase diagram of the GN model in the continuum from [Q| (left) and on the lattice 
(right). Full lines correspond to second order phase boundaries, dotted lines denote the (incorrect) first order 
phase boundary from the translationally invariant calculation. Pl is a multicritical Lifshitz point and C is 
the Euler constant. The dashed lines in the right figure illustrate the finite size artefacts, in particular the 
incommensurability effects at the right phase boundary of the crystal phase. 



or in momentum space 
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d 2 k 



(2.7) 



(2tt) 2 k 2 + a 2 

which is a self consistency equation for the fermion condensate or simply the gap equation. Equiv- 
alent equations can be derived via Hartree-Fock, Schwinger-Dyson or Bethe-Salpeter approaches. 
With the help of these equations one can derive the spectrum of the GN model ^ - it contains 
fermions with mass m = On> rc-fermion bound states, and baryons with mass m B = Oq ■ ™. 

What makes the GN model particularly interesting for our purpose is its rich phase structure 
in the (jU,r)-plane where pL represents the fermion chemical potential and T the temperature . In 
the homogeneous mean field approximation this was studied by Wolff [Q] who found two phases, a 
massive and a massless Fermi gas, separated by a line of first and second order transitions meeting 
in a tri-critical point, cf. Figure [j]. Recently the phase structure has been further clarified by Thies 
and Urlichs [Q, ||] who relaxed the tacit assumption of translational invariance of the condensate. 
Motivated by the fact that matter at low density forms isolated baryons they analytically found an 



inhomogeneous solution a(x) to eq.([23|). Indeed, they discovered that there exists a new baryonic 
matter phase at low temperature where baryons form a crystal structure. The transitions from the 
massive to the crystal phase and further to the massless phase are all second order. 



3. Lattice formulations of the Gross-Neveu model 



Let us now consider the GN model on a two-dimensional lattice, 

\2 N 



x ZA x,y a=\ 



where the staggered Dirac operator 

l r 

2 



5 x,y+l ~ 5 *,v-l 



x,y+2 x,y—2 



(3-1) 



(3.2) 



2 In the large-Af limit, a discrete symmetry can break spontaneously in (1 + \)d even at finite temperature. 
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describes 2 flavours and 



-8 xy (o(x) + a(x- 1) + a(x- 2) + a(x- 1 - 2)) 



(3.3) 



The modification of the naive discretisation a(x)8 xy — > L xy is necessary to ensure the correct con- 
tinuum limit [Qj. With the staggered discretisation the following discrete chiral symmetry is pre- 
served, 

X ( x )^(-iy^ x ( x ), 2 (x)^-(-iy^x(x), o(x)^-a(x). (3.4) 

Alternatively we consider a discretisation respecting exact chiral symmetry by employing the over- 
lap Dirac operator 



D = m\\+D w (-m) D^(-m)D w (-m) 



-1/2 



(3.5) 



satisfying the Ginsparg-Wilson relation + D = ^D^D. The coupling of the fermionic fields to 
the scalar field is introduced according to 



jgf = y(x) 



D x ,y 



Am 



D x ,y~D 



x,y 



4m 



o(x)8 xc 



¥(y) 



(3.6) 



which is consistent with a scalar density transforming covariantly under a lattice chiral symmetry 
transformation. For a homogeneous condensate a — > const, it is just the usual mass term of the 
overlap operator. 

The large-Af limit of the lattice theory is obtained by minimising the free energy. Using a 
homogeneous mean field, we obtain a as a function of X from the the gap equation. In Figure || 
we show these scaling functions for various lattice sizes. The dashed red curve describes the 
asymptotic scaling curve 2 ? 'l 2 e~ 7l l 2X for staggered and 1.5539...<?~ jr / A for overlap fermions. As 
expected, the staggered fermion formulation exhibits an additional factor of 2 in the exponent of 
the scaling curve due to the doubling of the number of flavours. 

To start our investigation of the ( i u,r)-phase diagram we first determine the chiral transition 
temperature T c at /I = where the chiral condensate a(r) vanishes. The results are shown in 
Figure |3| where we plot the scaling of T c /oq versus (aOo) 2 on various lattice sizes for the staggered 




Figure 2: Scaling of the condensate a as a function of the coupling A for staggered (left) and overlap 
fermions (right) for various lattice sizes. The dashed red lines are the asymptotic scaling curves. 
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Figure 3: Scaling of T C J Co vs (ado) 2 on various lattice sizes for the staggered (left) and the overlap operator 
(right). 

operator on the left and the overlap operator on the right. Both formulations exhibit lattice artefacts 
of similar size, but with a different sign, and nicely approach the analytically known value T c / Oq = 
e c j% = 0.5669... in the continuum [gp, thereby confirming the universality of the continuum limit. 

4. Baryonic matter in the Gross-Neveu model 

In both formulations a finite chemical potential can be introduced by weighting the temporal 
derivatives with factors exp(±ju) 0]. We can then check the universality of the continuum limit 
of \i c at T ~ and look at the discretisation artefacts for the two lattice formulations. In Figure 
|] we show the scaling of jU c /<7o versus (acb) 2 for the staggered Dirac operator on the left and the 
overlap on the right, still using only the homogeneous ansatz for the solution of the gap equation. 
Again, the two formulations show discretisation artefacts of similar size, but with a different sign, 



-- L=100 
-- L=200 
-- L=400 




Figure 4: Scaling of \l c j Co vs (aob) 2 on various lattice sizes for the staggered (left) and the overlap operator 
(right) using the homogeneous ansatz for the solution of the gap equation. 

and nicely scale to the analytically known value \i c = 1/ y/2 from the homogeneous ansatz of the 
condensate [Qj. However, in both cases the continuum value is approached non-monotonically - 
this is caused by the fact that the Fermi momentum changes continuously with jj. while the lattice 
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momentum is quantised by the finite lattice size. This finite volume effect would be most striking at 
exactly T = where one would expect a sawtooth behaviour; at T > 0, however, it is smoothened 
due to the softening of the Fermi surface. 

In Figure |5] we show the full phase diagram for the staggered operator at weak (L t = 80) and 
strong coupling (L t = 4), still using the homogeneous ansatz for the condensate. We note that 
below the tricritical point both the region of metastability and that of the chirally broken phase 
shrink considerably towards strong coupling. 
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Figure 5: Phase diagram at weak (L t — 80) and strong (L t = 4) coupling using the homogeneous ansatz 
for the condensate. The full line above the tricritical point marked by the dot denotes the second order 
transition, the dashed line below denotes the first order ttansitions while the two full lines to the left and 
right of it bound the region of metastability associated with the first order transition. 

Let us now relax the assumption of translation invariance. In order to determine the phase 
boundaries related to the transitions into the new baryonic matter phase, we perform on the one 
hand a brute force minimisation of the free energy and on the other hand check for instabilities 
of the homogeneous vacuum. The latter is achieved by monitoring the eigenvalues of the Hessian 
matrix in the homogeneous vacuum - a negative eigenvalue indicates that the free energy can be 
lowered further by an inhomogeneous perturbation. This is illustrated in Figure |6| where we show 
the lowest eigenvalue of the Hessian matrix associated with spatial variations of the condensate as 
a function of /I for a fixed temperature. It is clearly seen that for ji > 0.075, where the preferred 
homogeneous condensate is a = 0, an inhomogeneous one is favoured. In fact, a brute force 
free energy minimisation shows that this is also true for some range pL < 0.075, however finite- 
size effects cause a small free energy barrier between the two vacua. We also note that the non- 
monotonic behaviour of the lowest eigenvalue is an artefact of the non-commensurability of the 
spatial lattice size with the intrinsic length scale of the inhomogeneous condensate (compare the 
two volumes in Figure ^). In Figure |l| right we show the new phase diagram in the thermodynamic 
limit, at fixed L t = 80, to be compared with the continuum phase diagram on the left. The dashed 
lines show the effects from finite volume. In particular we point to the fringes at the right phase 
boundary of the crystal phase which are the reflection of the incommensurability effects discussed 
above. In the thermodynamic limit these effects disappear. 
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Figure 6: The lowest eigenvalue of the Hessian matrix associated with spatial variations of the condensate 
as a function of jj, for a fixed temperature and two volumes. 

5. Summary and Outlook 

We have investigated the phase structure of the GN model in the (/i, T)-plane. The breakdown 
of translational invariance of the chiral condensate requires a revision of the GN model phase 
diagram. Besides the massive and massless Fermi gas phase, a new phase of baryonic matter 
emerges which forms a baryon crystal. The transition to the new phase is always second order 
for any temperature. Our investigation of the phase structure on the lattice indicates that the new 
crystal phase disappears at strong coupling, because the topological excitations associated with the 
forming of the crystal fall through the lattice. Furthermore, large volumes are needed to detect the 
baryon crystal phase and to avoid or reduce complicated artefacts due to the incommensurability 
of the intrinsic length scale of the inhomogeneous condensate and the lattice volume, which distort 
both the phase boundary and the order of the phase transition. 

In this exactly solvable model, the crystal is formed by topological defects, the kinks and 
antikinks, which carry the baryon number. This crystal becomes stable for a sufficiently large 
chemical potential. One may wonder how general this phenomenon is. A similar topological 
crystal may occur in the (2+ \ )d Nambu-Jona-Lasinio model. If it does, perhaps it also occurs in 
QCD. 
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